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Lecture  5 
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STABILITY  OF  THE  PLAIIE  DEFLAGRATION  WAVE 

Steady,  plane  deflagration  was  introduced  in  the  second  lecture;  here 
ve  shall  consider  infinitesimal  perturbations  of  it  and  so  examine  its 
stability.  We  shall  find  two  basic  phenomena  -  the  hydrodynamic  and  Lewis- 
number  effects. 

Without  the  constant-density  approximation  our  task  is  not  easy,  because 
the  perturbation  equations  (though  of  course  linear)  have  variable 
coefficients.  There  are  three  ways  in  which  this  difficulty  can  be  overcome, 

(i)  If  attention  is  restricted  to  disturbances  whose  wavelength  is 
much  greater  than  the  thickness  of  the  deflagration,  a  hydrodynamic 
description  is  appropriate  (section  3.1)-  This  eliminates  the  Lewis- 
number  effect,  though  it  may  be  readmitted  as  a  perturbation  (Pelce 
&  Clavin  1982) . 

(ii)  If  the  con st ant -density  approximation  is  adopted,  then  the  T- 
and  '/-equations  (the  only  ones  that  have  to  be  solved)  have  constant 
coefficients.  This  eliminates  the  hydrodynamic  effect,  leaving  the 
Lewis-number  effect.  The  hydrodynamics  can  be  reintroduced  in  the 
context  of  a  weakly  nonlinear  theory  based  on  appropriately  small 
heat  release. 

(iii)  For  one-dimensional  disturbances,  the  governing  equations  can 
be  reduced  to  their  constant-density  form  by  means  of  a  von  Mises 
transformation,  the  distance  variable  being  replaced  by  the  particle 
function.  This  requires  the  diffusion  coefficients  to  be  proportional 
to  T  (rather  "char,  constants),  by  no  means  an  unreasonable  assumption 
physically. 

In  this  lecture  we  shall  consider  the  limiting  cases  (i)  and  (ii),  thereby 
isolating  the  hydrodynamic  and  Lewis-number  effects. 
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1.  Darrieus -Landau  Instability 

The  treatment  of  flame  stability  from  a  hydrodynamic  viewpoint  is  due, 
independently,  to  Darrieus  and  Landau.  The  work  of  Darrieus,  a  French 
aeronautical  engineer  well  known  for  his  invention  of  the  vertical-axis 
windmill  (Darrieus  rotor),  i3  often  cited  as  part  of  the  Proceedings  of  the 
19^*c  International  Congress  of  Applied  Mechanics,  but  these  were  never 
published.  Copies  of  a  1938  typescript  are  in  the  possession  of  several 
members  of  the  combustion  community. 

As  we  saw  in  section  3.1,  large-scale  disturbances  of  a  plane  flame 
are  described  by  Euler's  equations  on  either  side  of  a  temperature 
discontinuity,  namely 


=  0,  pDv/Dt  =  -Vp  -  pgi. 


(1) 


Pars  have  been  dropped,  so  it  should  be  remembered  that  the  scales  are 
much  bigger  than  the  diffusion  time  and  length.  Note  that  a  gravity  term 
has  been  added,  corresponding  (when  g  is  positive)  to  the  burnt  gas  in 
v  >  0  lying  above  the  fresh  mixture  in  x  <  0.  The  jump  conditions  (3.9,10) 
apply  across  the  discontinuity,  whose  nominal  position  is  x  =  0. 

The  undisturbed  flow,  found  as  a  solution  of  equations  (l)  satisfying 
the  jump  conditions  (3.9>10),  is 
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if  the  (small)  disturbance  parameter  is  absorbed  by  F^  (see  figure  l); 
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ve  consider  perturbations  for  vhich 

=  A  exp(iky  +  at)  with  k  >  0.  (4) 

Restriction  to  two  dimensions,  implied  by  this  and 

Y  =  (u,v) ,  (5) 

involves  no  loss  of  generality.  The  goal  is  to  determine  the  growth  ' 
parameter  a  as  a  function  of  k,  the  prescribed  wave  number  of  the 
disturbance.  The  corresponding  perturbations  of  the  flow  field  are 
governed  by 

-  0,  7*^  *  0,  p03y1/3t  +  p0u03y1/3x  =  -Vp^  (6) 

and  the  problem  is  complete  once  we  have  found  the  jump  conditions  satisfied 
by  these  perturbations. 

Since  only  terms  that  are  linear  in  disturbance  quantities  are  retained, 
the  -unit  vectors  in  the  directions  normal  and  tangent ial  to  the  discontinuity 


an  assumption  that  can  be  justified  for  SVFs  with  £  =  1  (section  U.2)  and 
for  ItSFs  (section  4.U).  The  conditions  (8),  (10)  may  all  be  applied  to  the 
undisturbed  location  of  the  discontinuity,  i.e.  x  =  0. 

Solution  of  the  perturbation  equations  proceeds  separately  on  the  two 
sides  of  the  discontinuity.  In  the  fresh  mixture  the  flow  is  irrotational 
(since  there  are  no  upstream  disturbances),  so  that 

t 

^  =  kP^e^,  v1  =  ikP^e^*  p1  =  -(a+k)Pfe“X  for  x  <  0,  •  (ll) 

if  the  factor  exp(iky+at)  is  omitted.  The  amplitude  P  of  this  potential 
field  is  as  yet  unknown.  In  the  burnt  mixture  there  are  also  rotational 
terms,  due  to  the  generation  of  vorticity  at  the  curved  flame,  so  that 

=  kP.^e  +kSe  ,  v^  =  -lxP^e  -(ia/a)Se  , 

p.  =  (a/o-k)P e  ‘  for  x  >  0.  (12) 

1  b 

In  addition  zo  P.  ,  the  amplitude  S  of  the  solenoidal  field  is  still  to 
b 

be  determined. 

The  requirements  (8), (10)  at  the  discontinuity  give  four  homogeneous 
equations  for  A,Pf,Pb,S.  These  have  a  nor.-trivial  solution  if  and  only 


(a+l)a  +  2aka  +  (a-l)(g-ak)k  =  0, 


-W-^k2  -  <S£)gk 


( a+l) 


In  the  absence  of  gravity,  i.e.  for  g  =  0,  the  larger  root  is  positive 
for  all  k,  corresponding  to  instability  for  all  wavelengths.  This  is 
the  Darrieus-Laniau  result,  which  has  been  an  embarrassment  to  combustion 
scientists  ever  since  it  was  discovered.  We  shall  now  see  how  gravity 


modifies  the  unacceptable  conclusion  that  plane  flames  are  unstable. 
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vfhen  g  r  0  it  is  instructive  to  examine  the  short  and  long  wavelength 
limits.  As  k  -*■ 


a 


l-a  ±  /a 


a{a  +a~ 1) ]k/(a+l)+. . . , 


(15) 


so  that  "he  influence  of  gravity  dies  out  and  short  wavelengths  are  unstable. 
As  k  -►  0 , 

a  -±t>/ (a- l)gk/(o+l)  -ak/(o+l)+. . .  (l6) 

so  that  gravity  stabilizes  long  wavelengths.  The  critical  wavenumber 
separating  stable  and  unstable  disturbances  is 

k  =  -Jo.  (17) 

Hydrodynamic  instability  is  observed  in  flames  (see,  e.g. ,  Sivashinsky 
1983),  but  it  is  often  absent .  There  are  several  ways  of  reconciling 
this  fact  with  the  present  results.  The  flame  may  be  too  small  to  be 
treated  as  a  hydrodynamic  discontinuity;  the  wavenumbers  that  allow  the 
disturbed  flame  to  be  so  treated  may  be  less  than  kc;  or  there  may  be 
other  stabilizing  influences,  such  as  curvature  (section  1). 

2.  Tne  lewis -Humber  Effect:  5V7s 

We  now  set  aside  the  hydrodynamics  and  investigate  the  Lewis-number 
effect,  by  adopting  the  constant-density  approximation.  The  first  pert 
of  our  discussion  is  concerned  with  377s,  which  are  governed  by  equation 
(1.12;  with  *?  a  0  if  volumetric  heat  losses  are  negligible. 

Section  3.3  found  that  such  flame  are  unstable  to  plane  disturbances, 
even  finite  ones,  for  L  >  1.  Tonplanar  disturbances  are  also  unstable, 
r.cv  for  L  <  1  as  well.  This  can  be  demonstrated  by  examining  large-scale 
disturbances ,  large  even  compared  to  -5.  To  this  end,  ve  introduce  new 


-5.6- 


(n,C,t)  =  6(n,C,r)  (18) 

and  consider  the  limit  5  -*■  0.  Then  3  =  0(5)  in  the  basic  equation 
(b.12),  and  W  =  1  to  leading  order;  to  next  order,  ve  find 

W  =  1  -jSbX  with  K  =  S-1dS/dT,  (19) 

an  explicit  formula  for  the  effect  of  weak  (superficial)  stretch  on  the 
flame  speed.  Clearly  such  stretch  decreases  the  wave  speed  only  for 
L  >  1  (cf.  section  k.3). 

The  instabilility  of  the  plane  flame  for  L  <  1  follows  immediately 
from  this  relation.  Suppose  that,  because  of  some  disturbance,  corrugations 
have  formed  (figure  2).  The  troughs,  as  viewed  from  the  burnt  mixture, 
experience  positive  stretch,  while  the  crests  are  compressed.  It  follows 
from  the  formula  (19)  with  b  <  0  that  the  flame  speed  at  the  troughs  is 
increased  that  at  the  crests  decreased,  so  that  the  amplitude  of  the 
corrugations  will  grow  (corresponding  to  instability.  Ve  are  dealing 
with  a  cellular  instability,  whose  nature  is  most  clearly  seen  for  the 
lTZ?s  discussed  below.  Because  of  the  connection  (3.27)  between  the  speed 
and  temperature  of  the  flame,  the  crests  are  relatively  cold  and,  hence, 
less  luminous  than  the  rest  of  the  flame.  For  L  >  1  the  effect  is  reversed, 
a  stable  situation. 

That  this  stability  conclusion  remains  valid  for  disturbances  of  more 
moderate  scale  can  also  be  seen  from  the  basic  equation  (1.12)  of  SVTs. 

Setting 

=  i+w  ,  :<  =  (20) 


>  (  m'..^ /d.x -<  L )  =  2'..\  with  W1  =  -F1t,  <1  =  F1 


nn 


(21) 
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x  =  -  t  +  ?1  (  n,  t) 


is  the  location  of  the  disturbed  flans  sheet.  Growth  of  the  mode  proportional 
to  exp(at-t-ikn)  is  therefore  determined  by 

2  — 1  ?  —1  /  —2  2* 

a  -2b  a-k  =  0,  i.e.  a  =  b  ±  /b  +k  .  (23) 

For  k  ^  0,  one  of  the  roots  is  positive  whatever  the  value  of  b  (i.e.  L) , 
so  that  there  is  instability  for  all  non-planar  disturbances  when  L  /  1 . 

For  k  =  0,  there  is  instability  only  when  b  is  positive,  i.e.  for  L  >1. 

The  instability  of  SVFs  restricts  their  value  as  a  framfcv  ~k  within 
which  to  discuss  flames ,  although  the  qualitative  insights  c  .ined  from 

such  solutions  can  be  very-  useful  (e.g.  concerning  flame  tip  see  Buckmaster 

f.'5J 

i  ludford  1982.V  In  pursuit  of  stable  flames,  the  more  usu?  lomenor. , 

*r 

ve  now  turn  to  NSFs. 

3.  The  Lewis-Number  Effect:  ?TEFs 

The  equations  governing  NEFs  were  derived  in  section  k.k.  To  obtain 
the  combusion  field  of  a  steady,  plane  deflagration,  we  introduce  the 
coordinate 


n  =  x  +  t  (2u) 

based  on  the  flame,  and  seek  a  solution  of  the  resulting  equations 

*?  7  ":  t+FT/3n-32T/3n2-32?/ay2  =  3h/ 3t+3h/ 3n-32h/3n2-^h/2y2-£(  32T/3n  +  32T/  ay)  =  0 

(25) 


depending  only  on  n.  This  yields 


W 


-  IX  „ne 


•  ho  = 


for  n  >  0 


ns  she  undisturbed  temperature  and  enthalpy  (perturbation)  profiles. 
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since  equilibrium  must  prevail  behind  the  flame. 

If  the  equation  of  the  disturbed  flame  sheet  (see  figure  l)  is 


n  =  ?1(y,t). 


(27) 


so  that  the  normal  derivative  in  the  jump  conditions  (k. 27-29)  becomes 
3/3n-?^y  3/3 y,  then  these  conditions  may  be  written 


*i  -  -Vr  fi(V  =  -nfFi*  3TI/3n  =  -WhiX 

6(3h1/3n)  =  ah^/lg-2*Yfr  with  tg  =  2T*/Yf, 


(28) 


all  quantities  being  evaluated  at  the  undisturbed  flame  sheet  n  =  0. 

(We  have  used  the  relation  =  -h1/T,o+...  and  superscripts  ±  to 
denote  values  at  n  =  0i.)  The  problem  is  to  solve  the  equations  (25) 
subject  to  these  conditions  and  the  requirement  that  T^,h^  die  out  as 
r.  -±* .  Arbitrary  initial  conditions  are,  as  usual,  taken  into  account  by 
considering  perturbations  for  which 


F1  =  A  exp(iky+at).  (29) 

The  possibility  of  stability  can  be  seen  most  easily  for  £  -  0,  when 
the  jump  conditions 


5(h)  =  6(3h1/3n)  =  0 


(30) 


ensure  that 


=  0  for  all  n  (31) 

is  she  appropriate  solution  of  the  differential  equation  (25b).  The  remainin 
problem  is  to  solve  the  T^-equation  for  n  <  0  alone,  subject  to  the 
boundary  conditions 
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T1  =  3^/an  =  -Y1A  at  r.  =  0; 


(32) 


factor  exp(iky+at)  has  bssn  omitted.  How,  the  solutions  of  the 
curbation  equation  are 

^(j—  $■-/-  <  =  /i+da+ltk2  ( — n /2  <  arg  <  £  t/2) 


(33) 


tere,  if  attention  is  restricted  to  unstable  modes  Re(a)  >  0,  we  have 


Re ( < )  -  1. 


(3»0 


follows  that  the  appropriate  solution  is 


Tx  =  S  e^1+l^n/2  for  n  <  0, 


(35) 


:i  then  the  boundary  condition  (32a)  requires 


3  =  *(!+<} 3,  i.e.  <  =  1  ora=  -k“ 


(36) 


:is  contradicts  our  assumption  that  the  node  is  unstable,  so  we  must 
■tcluie  that  there  are  no  unstable  nodes:  when  L  =  1  the  flame  is 
able  for  all  wavenumbers  k. 

7 ron  the  results  obtained  in  section  2  for  slowly  varying  disturbances, 
a  nay  expact  that,  as  becomes  large  (positively  or  negatively),  long 
veler.wth  disturbances  (k  small)  will  become  unstable.  To  see  how  this 


the  Jump  conditions  (28)  for  f,  /  0. 


lotion  o:  the  perturbation  equations  proceeds  separately  on  the  two 
des  of  the  flame  sheet.  In  front  we  find 

7.  =  re(l+<)r-/2,  =  {C+;<“1rk2-(l+<)2/U]Bn>e(l+’c)n/2  for  n  <  0,  (37) 

.ere  <  has  the  definition  (33b)  and  again  the  factor  exp(iky+at)  has 
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T1  =  °*  hi  =  De^1-*^2  for  n  >  0.  (38) 

The  expressions  are  valid  one -  for  Re(l+<)  >  0  and  Re(l-ic)  <  0,  i.e. 
when  she  condition  (3*0  is  satisfied,  as  it  is  for  unstable  nodes.  The 
jump  conditions  (23)  nov  yield  a  homogeneous  system  for  the  coefficients  ■ 
A,3,C,D;  a  non-trivial  solution  exists  only  if 

2tc2(l-<)  +  [(l-<)2-Uk232^  0  with  1=1/1,  (39) 

o  s 

a  result  due  to  Sivashinsky. 

This  dispersion  relation  should  be  viewed  as  determining,  for  each  £, 
the  growth  parameter  a  of  any  unstable  disturbance  mode  of  wave  number  k 
The  stability  boundaries  in  the  Z,  k -plane  therefore  correspond  to  Re(a) 
they  are  shown  in  figure  3.  At  any  wave  number  there  is  a  finite  band 
of  lewis  numbers,  always  including  L  =  1 ,  for  which  the  flame  is  stable 
and  outside  of  which  it  is  unstable.  The  left  boundary,  on  which  Im(a) 
vanishes  also,  is  a  creation  of  nor.planar  disturbances:  the  dispersion 
relation  (31)  does  not  yield  an  a  with  positive  real  part  for  £  <  -1 
and  =  0.  In  fact,  the  unstable  mode  becomes  neutral  as  this  part  of 
the  '.-axis  is  approached.  Or.  the  other  hand,  the  instability  predicted 
by  the  right  boundary ,  where  Tm.(a)  4  0,  does  survive  a  planar  treatment: 
the  dispersion  relation  (39)  with  k  =  0  does  yield  an  a  with  positive 
real  mart  for  Z  >  2(1+  /3). 

The  snabilility  boundary  in  the  neighborhood  of  the  point  Z  »  -1, 

>  =  C  can  be  explained  by  flame  stretch,  as  was  the  destabilizing  effect 
of  very  long  wavelength  disturbances  for  L  <  1  in  section  2.  We 
reei  rr.iy  r.coe  the  relation  (U.kg)  between  weak  stretch  and  flame  speed: 
for  i  <  -1,  positive  (negative)  stretch  decreases  (increases)  the  wave 
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speed,  tending  to  increase  the  amplitude  of  any  long-wavelength  disturbance 
(cf.  figure  2);  for  l  >  -1,  amplitudes  are  decreased. 

There  are  convincing  reasons  for  believing  that  the  left  boundary 
is  associated  with  cellular  flames,  a  common  laboratory  phenomenon: 
it  corresponds  to  L  <  1,  the  Lewis  numbers  for  which  such  flames  are 
seen;  the  crests  of  disturbed  flames  are  darker  than  the  troughs  (just 
as  for  SVFs  with  L  <  l),  this  being  a  characteristic  of  cellular  flames; 
and,  for  i  <  -1,  all  modes  with  k  <  W-il+l)  are  unstable,  suggesting 
that  the  outcome  of  the  instability  will  not  be  monochromatic,  smother 
characteristic  of  cellular  flames.  In  lecture  6  we  shall  give  a  nonlinear 
theory  that  arises  naturally  from  the  present  linear  analysis  and  reinforces 
this  conviction. 

The  right  stability  boundary  is  associated  with  pulsations  or  traveling 
waves;  it  is  relatively  inaccessible  because  L  is  rarely  much  bigger 
than  1.  However,  similar  phenomena  may  be  expected  in  so-called  thermites, 
for  which  L  =  °°.  Otherwise,  special  means  must  be  devised  to  make  the 
boundary  more  accessible.  These  matters  form  the  subject  of  lecture  T. 

h.  The  Role  of  Curvature 

The  olane  premixed  flame  whose  stability  has  been  discussed  so  far 
is  an  idealization  seldom  approximated,  since  in  practice  the  flame  is 
usually  curved .  Even  under  -circumstances  designed  to  nurture  a  plane 
state,  imperfections  can  thwart  the  best  efforts  and  give  a  curvature, 
albeit  weak,  to  the  flame.  In  this  section  we  shall  investigate  certain 
slightly  curved  flexes,  those  amenable  to  the  SVr  analysis  in  section  2 
and  those  associated  with  the  left  stability  boundary  for  NEFs  in  section  3. 

Consider  the  cylindrical  source/sink  i'Iow 

u  =  iR/r ,  v  =  C,  (-0) 
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where  9r  is  radial  distance  and  u,v  are  now  polar  components.  The 
undisturbed  flame  is  then  circular,  with  9-multiplied  curvature 


K 


(1*1) 


determined  by  the  requirement  that  the  mass  flux  through  the  unstretched 
flame  is  1.  The  expressions  (20)  are  replaced  by 

W  =  1  +  W1,  K  =  ^  +  kW2,  ‘  (U2) 

where  is  due  to  the  displacement  of  the  flame  in  the  source/sink 

flow,  which  (for  the  constant-density  approximation)  is  undisturbed. 
Substitution  in  the  basic  equation  ( 1* .  12 )  gives 


b ( dW  /dx  — V.v  —  K, —  kW  )=  2W 
1  "_L~  11-  1  1  1 


(1*3) 


with 


x  =  -(R"1F1±F1X) ,  2  -X, 


■+  2R"2?1  , 


<1  =  ±r~2(f1+f1(j^)  (1*1*) 


r  =  R  +  F1(<j>,r)  (**5) 

is  the  location  of  the  disturbed  flame  sheet  (figure  l).  The  behavior  of 
the  mode  oroportional  to  exp(ctt  +  in$)  is  therefore  given  by 

a2-2(b"1;R1)a-R“1(±2b"1+n2R_1)  =0.  (1*6) 


To  discuss  this  result  we  introduce  the  wavenumber 

k  =  R_1n  (1*7) 

and  note  that  the  equation  (23a)  is  recovered  as  R  -*•  <*>.  There  will  be 
instability  for  those  values  of  k,  R  that  make  one  of  the  (always  real) 
roots  of  ecuation  ( Uo)  positive.  Such  values  are  bounded  by  ones  for 
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which  a  =  0,  but  the  reverse  is  not  true.  Figure  4  shows  the  line  on  which 
one  root  vanishes  and  the  corresponding  value  of  the  other  root;  the  full 
part  of  the  line  is  therefore  a  stability  boundary  but  the  dashed  part 
is  not.  Results  for  ?-*■«,  given  in  the  discussion  of  equation  (23), 
then  shov  that  the  stability  regions  are  as  lableled.  We  conclude  that 
curvature  is  never  stabilizing  for  a  sink  flow,  but  that  it  is  for  a  source 
flou  when 


L  <  1,  k  <  /-2/bR. 


We  turn  nov  to  the  neighborhood  of  l  =  -l,k  =  0  in  figure  3,  where 
the  dispersion  relation  (39)  reduces  to 

a  =  -(l+7)k2-4kU;  (1) 


:e  have  a  balanced  equation  ve  must  require 


k2  *  1(1+*.),  so  that  a  =  0(l+i)2. 


To  see  how  this  relation  is  modified  by  curvature,  ve  consider  once 
more  a  circular  flame  r  =  R  sustained  by  a  source/sink  flow.  A 
derivation  of  the  new  dispersion  relation  will  be  given  in  section  6.2, 
where  it  is  needed  for  a  nonlinear  analysis.  Here  ve  shall  only  note 
that  the  connection  between  wave  speed  and  stretch,  represented  by  the 

generalization  (!?'  of  the  relation  (4.U8),  is  preserved;  but  that 

2 

recalculations  of  the  stretch,  represented  by  the  term  in  k  ,  and  of 
the  wave  speed,  represented  by  the  term  a  ,  are  needed  because  of  the 
non-uniformity  of  the  velocity  field. 


The  modified  dispersion  relation  is 

a  *  -(l+Tjk2-!*1*  ;  R~\  (50) 

vhere  k  is  now  limited  to  the  values  (1*7);  a  balanced  equation  requires 

R  =  0(1+I)~2  (54J 

in  addition  to  the  earlier  restriction  (50a).  For  Z  <  -1  and  no  curvature, 
there  is  a  band  of  wavenumbers 

0  <  k  <  }/-(l+z)  (52) 

for  which  a  is  positive,  corresponding  to  instability  (cf.  figure  3). 
Source  flow,  i.e.  convexity  of  the  flame  sheet  towards  the  burnt  mixture, 
narrows  the  band  (figure  5)  and,  for  sufficiently  high  curvature 

R"1  >  (l+?)2/l6,  (53) 

eliminates  if.  Sink  flow,  on  the  other  hand,  widens  the  band  indefinitely. 
For  i  >  -1  and  no  curvature,  a  is  negative  for  all  k,  corresponding  to 
stability.  Source  flow  does  not  change  this  conclusion,  but  sink  flow 
produces  instability  for  wavenumbers  below  an  ever-increasing  upper  limit 
(figure  5) . 

For  both  the  SV?  and  the  NEF,  convexity  towards  the  burnt  mixture  is 
found  to  be  stabilising,  concavity  destabilizing. 
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Figure  Captions 


5.1  Notation  for  stability  analysis  of  plane  flame  as 
hydrodynamic  discontinuity  (x)  and  as  NEF(n)  . 
Flame  surface  is  given  by  F^(y,t)  . 

5.2  Instability  of  SVF  for  L  <  1  . 

5.3  Linear  stability  regions  for  plane  NEFs.  The 
boundaries  are  P:  4k2  ■*  -(1+1);  B:  X  -  2(l+8k2) 
[1  +  /3(l+8k2j]/l+12k2). 

5.4  Effect  of  curvature  on  SVF  stability. 

5.5  Effect  of  curvature  on  NEF  stability. 
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